


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 


DSpace Repository 


Theses and Dissertations 


1974 


1. Thesis and Dissertation Collection, all items 


Heave dynamics of a surface effect ship. 


Dewey, William Alonzo. 


Massachusetts Institute of Technology 


http://ndl.handle.net/10945/17047 


Downloaded from NPS Archive: Calhoun 


ai | pune 
mM \ LIBRARY 


http://www.nps.edu/library 





Calhoun is the Naval Postgraduate School's public access digital repository for 

research materials and institutional publications created by the NPS community. 

Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
appointed — and published — scholarly author. 


Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 








~ 


HEAVE DYNAMICS OF A 
SURFACE EFFECT SHIP 
by 
WILLIAM ALONZO DEWEY 


Bevo. University of Utah 
(1968) 


SUBMITTED IN PARTIAL FULFILLMENT 
OF THE REQUIREMENTS FOR THE 
ENGINEERS DEGREE IN OCEAN ENGINEERING 
AND MASTER OF SCIENCE IN MECHANICAL ENGINEERING 
at the 
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 


June, 1974 





Minas KNOX LIBRARY 
“OSTGRADUATE SCHOOE 
LY, LtFC oma °39aq 


fv 


HEAVE DYNAMICS OF A 


SURE ACE EPR ECT SHIP 


by 
William Alonzo Dewey 


Submitted to the Department of Ocean Engineering and 
Department of Mechanical Engineering on May 10, 1974, in 
partial fulfillment of the requirements for the Engineer 
Degree and Master of Science. 


ABSTRACT 


The dynamics of the sidewall-type surface effect ship 
(SES) for vertical plane heave motion in response to waves 
are established through theoretical model and experimental 
investigation. The theoretical models are digitally simu- 
lated and compared to an experimental tow tank model SES 
(length to beam ratio of 6.5) for verification. The exam- 
ination of transient and frequency responses provides 
insight into the fundamental dynamic characteristics. 
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NOMENCLATURE 


AL Waterplane area of cushion, £-°b. 
A. Exit aimilow ames, 
A ap A os Bow and stern exit air flow areas, 
A Entrance air flow area. 
a Sidewall waterplane area, 4-y... 
a Added mass, 
B ieain Overall, . 
b Cushion width. 
b. Cushion damping coefficient. 
ibe Sidewall damping coefficient. 
b, db. + Digg 
Co: C. Quadratic orifice coefficients. 
BF External applied farce. 
g Gravity constant. 
h Exit cee melclt, 
ho Equivalent cushion height at equilibrium. 
h, Equilibrium gap height. 
k Polytropic constant. 
kK Sidewall stiffness coefficient. 
Ko Equivalent cushion stiffness coefficient. 
kK, Equilibrium buoyancy force on immersed sidewall, 2p Meu 
is Length of SES overall. 
L Length of cushion. 
*M, Air mass in cushion region. 
*M Supply flow rate. 


Aa equilibrium conditions are indicated by a zero subscript. 
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NOMENCLATURE (continued) 


Air mass flow rate into cushion. 

Air mass flow rate out of cushion. 
Mass of the craft. 

Ambient air pressure, absolute. 
Cushion air pressure, absolute. 
Supply air pressure from fan, absolute. 
ae. 

Craft forward velocity. 

Cushion volume. 

Duct volume. 

Volume of immersed sidewall at equilibrium. 
Wave frequency. 

Encounter frequency. 

Natural frequency. 

Horizontal position along craft length. 
sidewall thickness. 

Wave height at any point. 

Weriieal height of crait. 

Cushion height. 

Water height inside cushion at bow. 
Water height inside cushion at stern. 
Sidewall height. 

Average wave height inside cushion. 
Damping ratio. 


Wave height. 


“AU equilibrium conditions are indicated by a zero subscript. 





NOMENCLATURE (continued) 


Wave length. 

Water density. 

Air density. 
Cushion air density. 


Supply air density. 
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CHARA Ro 1 
INTRODUCTION 


Pea Backcround 

In the past few years substantial interest has developed in surface 
effect ships (SES), specialized waterborne air cushion vehicles which 
can provide high speed waterborne transportation. The SES configura- 
tion being studied is of the captured air bubble type. It employs the 
air cuShion vehicle principles, utilizing a pressurized compartment of 
trapped air to support the majority of the craft's weight. A smaller 
fraction of the total support force comes from the hydrostatic, hydro- 
dynamic and aerodynamic forces. The pressurized compartment of air 
(cushion) is contained by rigid sidewalls along the craft length and by 
seals both fore and aft. See Fig. 1.1. 

In the ocean environment, the principal vertical plane motions of 
interest from the aspect of seakeeping and habitability are heave, surge, 
and pitch. So as to provide reasonable estimates of these motions, 
effort is directed towards the development of tools that will provide 
information on the expected behavior of the full-scale craft. The 
principal tools uSed here are experimental model tests and mathematical 
Simulation via digital computer, 

The dynamic simulations are helpful where known Froude scaling 
laws do not hold. Since this is true of the heave-pressure coupled 
motion for the SES, then the simulation model becomes a very useful tool 
for full-scale motion prediction. This does not preclude the use of 
experimental scale models but, rather, it means that once digital 
simulations accurately predict experimental model results, then this new 


tool can be easily adjusted for known parameter changes. Therefore, it 
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can predict more accurately full-scale craft dynamics. 


1,2 Scope of Thesis 


The principal thrust of this thesis initially leads in two directions. 
These are (1) to develop a mathematical model of the equations of motion 
for digital simulation and (2) to build and test a scale model of an SES in 
the MIT tow tank facility. Attention is focused to the single degree of 
freedom in heave only. Heave motion of the SES, being the main contrib- 
utor to the craft vertical acceleration levels, is a good Starting point for 
digital simulation and experimental verification. 

Chapter II explains the modeling of the fluid suspension system. 
The primary components to the model are developed for later use in the 
Simulation Models A, B, and C found in Chapters III and IV where: 

Model A is the SES craft transient performance analysis on 
cushion with no entrance or exit air flow. (Principally this is an 
inverted open-bottom container suspended by a pressurized air mass on 
the water surface, ) 

Model B is the same as Model A except, now, air flow is considered, 

Model C is the simulated responses to regular waves, 

The most basic model consists primarily of the following major 
components: (1) a cushion region formed under the craft that is main- 
tained by low pressure air for craft support; (2) a nonflexible 
sidewall structure; (3) the forward and after seals; and (4) a fluid 
source for feeding air to the cushion. : 

The simulated transients and wave frequency analysis are to be 
compared to experimental model tests for verification, Comparisons of 
heave amplitude, phase and acceleration levels are illustrated, along 


with their time histories at selected frequencies, The experimental 





model tests were conducted at zero forward speed. 


1.3 Prior Work Related to Heave Motions 

Historically there has been interest in the analysis of air suspen- 
sion systems in heave from two major technologies: (a) the land-based 
air cushion vehicles, and (b) the sea-based air cushion craft. It is the 
latter field to which this thesis addresses itself, although useful informa- 
tion and analysis has been obtained from both areas, 

The most recent and extenSive analysis and simulation studies for 
SES craft have been conducted by Oceanics [1 Ie Simulation model 
analysis of the SES vertical plane dynamics has been performed, though 
verification is not extensive in the report. 

Other land-based air cushion suspension studies have been 
conducted for the Department of Transportation at MIT [2 and 3], the 


results of which have been widely published. 





ok 
[ ] indicate references. 


Note: From henceforward, "model" and 'modeling" shall refer to the 
theoretical model for simulation. In order to avoid confusion in 
terminology, ‘experimental model" or "SES experimental craft" 
is used when referring to the tow tank test craft. 
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Charli i 


THE HEAVE SUSPENSION MODEL 


2.1 The General Suspension Configuration 

A nonlinear dynamic model for heave motion of a rigid body is 
developed for the configuration illustrated in Fig. 1.1. The basic sus- 
pension system consists of (1) cushion region, (2) supply system, 
(3) seal system, and (4) sidewalls. Inthe configuration a fan supplies 
flow through a constant volume duct ane aves restriction to the cushion 
chamber (plenum). The flow from the plenum escapes through the 
Sealing gap area formed between the seal and the water surface. The 
fluid suspension is basically a three energy port device where energy 


transfer between the waves, air Suspension and craft takes place as 


illustrated in Fig. 2.1. In developing models for these parts of the 


Gravity 


Heave 
Craft 
Acceleration 
Forces 
and 
Motions 
Feed ; Leakage 
System —™ SUSIDQOS ee System 


Forces 
and 


Motions 





Water Waves 


Fig. 2.1 Craft - Suspension - Wave System. 





system, the following assumptions are made: 

(1) The cushion and source absolute pressures are comparable 
to the absolute ambient pressure. For the experimental 
tow tank model the cushion pressure is on the order of one- 
tenth psig and the supply pressure of one to two psig. 

(2) The external flow and pressure fields induced around the 
vehicle while it 1s moving with a forward velocity are not 
considered to affect the exit flow from the cushion. This 
assumption is not valid at high speeds, especially at the 
bow seal, Compensation for this effect could be modeled 
into the system for more extensive study. 

(3) The cushion chamber base area Ay is constant. 

(4) The pressure P. and density p. in the cushion region are 
Spatially uniform. 

(5) The Froude-Krylov hypothesis assumes the wave system 
iS not appreciably altered by the presence of the craft. 
Thus, it is assumed that the wave is not diffracted or 
attenuated as it passes under the cushion chamber. 

(6) The craft is assumed to be a rigid body. 

The principal forces acting are identified in the rigid body 


equations of motion for heave, Fig. 2.2. 


asl 
(m+a) | = 4+ 4, eep (1) 
where 
Ze tOnees Gue tO the cushion refion 
‘ (la) 


I 


“25, (P.-P,) ween 


16 








(a) Coordinate System. 


Craft Weight 


Hydrostatic and Dynamic Wave Forces 
(b) Cushion Forces, 
Craft Weight 





Hydrostatic and Dynamic 
| Sidewall Forces 


(c) Sidewall and Craft Forces 


Stern Seal 


Pao Cs c 


Hydrostatic and 
Dynamic Forces 
Bow Seal 


(d) Bow and Stern Seal Forces. 


Coordinates and Principal Heave Forces. 


iv 








where 


SW 


sidewall forces 


. ee (1b) 
“1S NA a 
forces due to Seal (considered small (1c) 
compared to other external forces) 
gravity and externally applied force 

(1d) 


mg +R, 


added mass 

waterplane area of cushion base 
waterplane area of the sidewall 
cushion damping coefficient 
sidewall damping coefficient 
external force acting on the SES 
stiffness coefficient of sidewall 
ie Acw 
Verelcalsmuovant Moree of iminersed 
Sidewall = 29 ve 

mass of craft 

ambient pressure absolute 

cushion pressure absolute 

craft position 

average external water elevation 


volume of immersed sidewall at 


equilibrium 
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Therefore, substituting (1, a-d) into a 
dz, 0 BS 
re rae a pb Oenes) | O77 > Mi @ps ee yt ME * FY (2) 


where 
b, 5 beet De 


In the analysis and later in the experiment, the forces acting on 


the seals were assumed to be relatively small compared to the cushion 


forces, mg + Er 


1 





Pa 
allt 
External . 
Water Internal \ 
Surface Water aa 
Cushion 
Surface Resi 
egion 


2.2 The Cushion Region Model 


Within the plenum chamber, conservation of mass is required; 
thus, the inlet supply flow (M,,) less the flow through the exit gap (Mot) 


must equal the rate of change of mass within the chamber: 


dM . ; ; 
dt . May : Nae (3) 
M, = air cushion mass 


The pressure in the plenum chamber may be related to the mass 
(M,) and chamber volume Sue) by the equation of state. A polytropic 
relationship is assumed for the cushion chamber. 


“ () indicate equations. 


ioe 
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p = (M./V)* constant (4) 


The polytropic process is felt to be between isothermal (k = 1.0) and 
isentropic (k = 1. 4). 

The total cushion volume ie at any instant in time may be 
expressed as the product of the base area Ay and the sum of the plenum 


equilibrium height es craft elevation Zy: and the water surface elevation 


Zo: 
Ue eee) (5) 
where 
ho = equivalent cushion height at equilibrium 
Zz. += average wave height under cushion 
5 Static depression of water plus the external 


wave elevation = (p.-p.)/P co 


2.3 Cushion Flow Entrance and Exit 

From (3) the entrance flow M.. and exit flow ieee are modeled as 
follows. First, for purposes of the experimental craft verification the 
inlet mass flow M,_, is considered constant since the entrance flow came 
from a high pressure source where the dynamics of the cushion region 
did not affect the mass flow to any degree. Secondly, in modeling an 
actual system where fan and duct dynamics are important, the mass flow 
into the cushion is modeled according to the classical quadratic orifice 
equation for non-compressible flow. Non-compressible flow may be 
assumed for the absolute source and cushion pressures Since they are 
comparable to the absolute ambient pressure (eee oem <1 psi). 


Mass flow into the cushion is: 
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. al 
Min = CoAgl2(ps-PQ) Py /p)” (6) 


The mass flow that exits from the plenum must either reverse and 
pass through the entrance orifice in the case where on @ Pe or the flow 
passes under the bow and aft seals. Again the Mout is modeled 


Similarly to the entrance mass flow with the exception that now the exit 


leakage area A, is varying. 


ay 
2 


Mout cael be Pa) ee /_? (7) 
Piere 
Aa = entrance flow area to cushion 
cn variable exit flow area under 
seals 
Co: C. = discharge coefficients 

(assume . 0.6) 

ye = fan supply pressure 


2.4 Exit Air Flow Region 


The exit area is the opening between the seal and the water surface 
through which the relatively low pressure air eScapes at the bow and 
Stern. This assumes that the sidewalls allow no leakage. The total 


leakage area is 


eu Ne | (8) 
e eb es 
where 
A. = total exit area 
Png LN = bow and stern leakage areas 
eb es 


The bow and stern exit areas are related geometrically to the craft 





seals and water vertical position by: 


b(-z, aera 1) if >0 
A. = (9) 
, 0 if <0 
L. 
= ‘f > 
b( Zt 2 at he) if >0 
Ans = (10) 
0 if <0 


where the wave height at the bow and stern respectively are given by: 


N 
! 


cushion water height at bow 


ob 
ae cushionmy ater height at stemn 
hy = equilibrium gap height 


2.5 The Fluid Source and Feed System 


This section on the fluid source is presented for future detailed 
Studies where the duct and fan dynamics may influence the full-scale 
craft results. The fan and duct dynamics, as given here, were not 
needed for the verification of the experimental model tests. The use of 
a constant mass flow rate source in the experimental test apparatus made 
it feasible to focus attention directly to the fluid suspension (cushion) 
System. The feeding system which supplies air flow to the cushion in 
the full-scale craft consists of a fan and a feed duct. The fan may be, 
in general, assumed to operate at a constant supply pressure. This 
assumption is valid if the source pressure remains nearly equal to the 
fan operating point equilibrium pressure. However, for large deviations 


from the operating point, the change in supply pressure as a function of 


ae 





Supply flow should be considered. 
The fan curves identified in Fig. 2.3 may be represented by the 


nonlinear functional relation: 


M,. = fp ,-P,) (11) 


Operation about the equilibrium point may be linearized for simulation 


and analysis wherein the mass flow rate M, is given as; 


M 


P. ~ P 
= = -5|/5— 4} +2 (12) 
where 
M, = supply flow rate 
M, = equilibrium supply flow rate 
Pee exelbins Norouhl aa, Shaloyophks she Sstabh ce 


The supply feeding duct is assumed to act as a fluid capacitance 


with duct resistance lumped into the plenum inlet supply area resistance. 


The equation of state for a polytropic pressure density relation- 
Ship similar to (4) combined with the conservation of mass may be 


written as: 





Vv 
a Ps = Mp- Mi, see 
Where 
Pp, = Supply air density 
Me = duct volume 


2.6 External Disturbances 
The transient performance study for Models A and B considers an 


instantaneous step increase in force F, applied to the model craft. In 


1 


7 
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the later study of Model C the external disturbance comes from the 
waves Zz. The external wave disturbance acts to change the air volume 
and hydrostatic force on the sidewalls, 

The typical seaway considered in the analysis is a regular Sinus- 
Oidal wave of height ie and frequency WwW. Inthe SES coordinate system, 
the origin is considered at the longitudinal center of gravity. Therefore 


the regular wave elevation for head Seas is described by: 








Ze = al. sin (<7 + wt) (14) 
where 
£ = £Jlength between Seals 
ee = wave height 
h = wavelength of wave 
u = craft forward velocity 
WwW. = encounter frequency = wt + we u/g 
for head seas 
W = wave frequency 
x = position along craft length 


from midship positive forward 
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The instantaneous wave height along the length of the cushion and 
at the seals is important in the accurate model representation of the real 
System. Inthe first case the water height at the bow and stern seals 
actS as an external disturbance that changes the equilibrium gap height 
he caused by air flowing from the cushion under the seals. The wave 
height at the bow and stern seals, without taking into account the 


depression due to cushion pressure, is given by: 


N 
1 


No sin (wt + 7l/) (15) 


N 
| 


n,Sin(wt - mL/d) (16) 


The average wave height under the cushion is determined by 
integration along the cushion length. The average wave height under the 


cushion is given by: 


n, £/2 
Zee = \ sin(2 7x/+ wt) dx 
: ~2/2 
Therefore 
n, 
ee ae sin(72/) sin w ot (17) 


This model is representative as long as the air pressure inside the 
cushion is considered a constant throughout and there is no localized 


wave pumping. 


2.7 Combined Model Equations 

The system equations defining the heave dynamics of the craft for 
Models A, B, and C in Chapters III and IV are given below, This 
includes the equations from sections 2.1 - 2.5. 


1 
(m+a)dz, /dt = -b,2, a8 Zee (pep) ak me + 


(2) dz, /dt = Z 





24 
(3) dM /dt = M,,- M 
(4) p. = (M._/V 8 constant 
Cc cfc 


(5) ve = AL(- Z4 + ho oa Zz) 


where Zo (p. - p,/ptz. 


° ad 
(6) M. = CA (2p. - Po)P5 fg)” 


. i 
a) Mout cee ee ; Pa)Po fx? 


(18) 
(8) A, = Alay aaa 
(9) Aap = b(-z, + Zot + ~“ it -0 
S 0 bi 0) 
(10) Aas = b(-z, a Pe h,) if >0 
0 } if <0 


(11) My = £(p, - p,) 


we (V 4 P5/KP 59) Zone Mr - Mn 


2.8 Linearized Model Equations 

A linearized set of equations for the heave mode was developed and 
used along with the nonlinear equations (18) for determining effects of 
cushion parameters. The linearized equations are given in state 
variable form for Models A, B, andC. Excluded from these three 
models are the duct and fan dynamics Since the experimental SES was 


Supplied from a constant mass flow source, 





eal 


where 


where 
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0-0 0 Bee ee ee eae 
=_— — Wen NSO Bo NS Nw 
Q t= (Wao wo == wo dO = 


a" 
bo 


= ~ARG 


1 0 A zy 0 0 
AF, 
Aga Ag3! {44,1 + | Bo, Boo 
TAA 
W 
0 | Az. A M. 0 Byo 
(19) 
A Z4 
: AF 
0 “18 A Z4 oa 0 Di9 
A Z 


= -(K, + K,)/(mta) 
= ~b, /(m+a) 


pC13/(mta) 


= ~K3C), ~ Kyl +7770) 


= -K3C,3 ~ K4C,3/p 


= 1/(m+a) 


“A 


- (k, - A,D,.)/(mta) 


] 
= -K,D,, - K,(1 + D,9/p) 


See 
= dp./dz, = p/(hoe/kp,, + 1) 
= op. /oM, SVB JO NUL iy (Foy iss) a9 38) 
= dp. /oz, = -(p/(hye/kp.g t+ 1) 


1 
2 


J 


on 
Cae eip! A,/l2(Po9 72) 
5h 
Co Pe/g) 2b(p. > Pp.) 


An@ 14 
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SURFACE EFFECT SHIP TRANSIENT RESPONSE 


3,1 Introduction 

The basic equations, as developed in Chapter 2 and summarized in 
2.7 and 2,8, are considered in the study of the transient response to 
external inputs. The theoretical results of various responses are com- 
pared to experimental data for verification. Two theoretical models 
tnat closely repeat the experimental eoststare designated by the titles, 
Model A (without air flow), and Model B (with air flow), As aresult of 
the transient responses, the natural frequency and damping of the SES 
are obtained, 


Prior to reading in detail the analysis in Chapters 3 and 4, turn to 


Appendix A for discussion of the experimental model and testing. 


3.2 Transient Performance Parameters 

In the transient performance analysis, linearized cushion model 
performance 1S compared to experimental test data. Response 
Pier actenictics are Cetermined tor several types of disturbances which 
are explained in 2,6. 

The Model A equations are representative of the heave dynamics of 
the SES on the air cushion without air flow passing through the boundaries 
of the craft. Basically, the Model A dynamic equations are Second 


order wherein the craft vertical position z, and velocity Z, are the 


1 
primary State variables. The cushion pressure and other desired out- 
puts are dependent variables of this system. The results of this 


analysis are the determination of the craft natural frequency and damping. 


The experimental data transients are conducted while the craft remains 
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stationary in the horizontal plane. 


From (18), Model A nonlinearized equations are given with the 


following conditions: 


dM ,/dt, dM,/dt, dp _/dt, Mee We = 


Upon linearization, the resultant equations in state variable form 


are listed in (19) where: 


Agog 431 Agg: Bgg, Cig = 0 


The second order system's characteristic equation in linear form 
is: 


(mt+a)Z, + b,z, + (K, + Ay e/(n, p/kp., + 1))z, = 0 (20) 


The quadratic equation indicates a natural frequency for the coupled 


heave-pressure mode given by: 


~k, + A, /(h_/kp__ +1) 
oe j —+-—_2 2 90 (21) 
n mta 
The damping ratio (relative to the critical damping) is given by 
— = b, /2w, (m-+a)) (22) 


In (20) the numerator terms are a measure of the system's stiff- 
ness, The first term kK, is the stiffness coefficient of the sidewalls. 
The second group of terms of the numerator is the stiffness coefficient 
due to the air cushion, It should be noted that the theoretical stiffness 
coeificient is approximately equal to the waterplane area of the cushion 
and sidewalls times the water density. 


All parameters of the second order system can be determined 


directly (See Table A,1) except for the damping b, and added mass a, 
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Estimates of these two variables may be derived from comparison of the 
transient heave response of the model and experiment. AS a check on 
these results from theory [5], estimates of the damping and added mass 
for oscillating pressure fields on a free surface, as formulated by 
stoker (1957), were used. Curves for added mass and damping from 
[5]are found in Appendix B. Note that, in order for correlation between 
experiment and theory to be valid, the seal and sidewall effects must 
also be estimated since 19% of the total’ waterplane area is due to the 
Sidewall. 

Figure 3.1 is the response for a step input in force acting on the 
craft. This force was 37.5 lb, equal to 100% of the craft weight. The 
natural frequency and critical damping ratios, calculated by comparing 


the digital response and experimental data, are approximately: 


® = 9,2 rad/sec 
mean 0 ae 
where m = 1.16, a = 4.34 (1b sec“ / ft), Ky + Ko = 461 15/ ft, b, = 30 lbsec/ft. 


Determination of the added mass, a, from theory based upon a 
two-dimensional oscillating pressure field [5] was in agreement with the 
experimental results. Theory from [5](see Appendix B), indicates 
that the damping ratio due to the cushion is on the order of €= 0.3. 

This appears to agree well enough with the experimental results, though 
some shecreoemes due to the sidewalls is certain to be present. 

Reference [1] indicated, though, that the damping due to the pressure 
field is negligible, which is contrary to what [5] and experimental results 
indicated. 


In this study it was expected that there would be more of a 
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difference between theory, based solely upon the cushion dynamics [5], 
and the experimental results, This is because the experimental craft 
had solid sidewalls and seals which covered 19% of the total waterplane 
area affected. But, within the accuracy of the experimental results, no 
detectable difference was found. Some theoretical estimates of the 
hydrodynamic force terms for the sidewalls could be found through the 
use of slender body theory, but for this thesis no Such estimates are 
made. Future studies should look into this. 

ippiecal velocity, acceleration and pressure Reanerent Weoeross 


ponding to Model A are shown in Fig. 3.2. 


oS. 0 Model B Transient Response 


The transient reSponse analysis for Model B is performed in the 
Same manner as in Model A. The principal addition to the dynamic 
equations is the entrance and exit mass flows described in (3). The 
entrance mass flow is assumed to be a constant while the exit mass flow 
is described by the quadratic orifice equation (7) where the exit area and 
cushion pressure vary with time, 

Again the resultant natural frequency Ww, and damping ratio € are 
determined, The state variables of the third order system are heave Z,, 


velocity z,, and air cushion mass M,. The resultant output variables 


le 
that are observed are the cushion pressure P. and acceleration Zs. AD GNE 


model equations from (18) have the following conditions imposed: 
dM ,/dt, dp. /dt = 0 and M,, = const. 


Upon linearization, the resultant equations in state variable form 
are listed in (19). These equations are used in the transient response 


study. 
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Prior to comparing model Simulations to test data, some observa- 


tions of the experimental tests are made, 


3.4 Transient Response from Experimental Tests 

The SES experimental model response with air flow through the 
seals presented some interesting observations, Figure 3.3 shows two 
transient response curves for two test conditions nearly identical except 
forthe seal design. Since the seals used in the experiments were non- 
flexible, the mass air flow rate would be expected to vary more dramat- 
ically here than for a deSign using a flexible seal that contours the water 
Surface, 

In test 1 of Fig. 3.3, the heave response shows a flattening at the 
top of the curve when the craft rises to 0.02' above equilibrium. In 
test 2 the craft rises to 0,038' above equilibrium. Why the difference? 

In test 1 the seal design was such that a vertical gate seal of $"' 
thick plexiglass was used to contain the cushion at the bow and aft. In 


this case the exit gap length was very short, resulting in a more rapid 


bleed-off of air. Intest 2 the seal was much thicker (2 inches). This 


Me 





Seal design Seal design 
forest 1. for test 2. 
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seal is faired to allow the water to flow smoothly around it without 
creating excess turbulence. The air escaping under the seal has to 
travel a larger distance before escaping from the plenum. For this 
reason the air loss is less, and the mean pressure level inside the 
cushion is higher. 

As a result of the seal design, the heave response of test 2 has a 
higher overshoot than in test 1. The corresponding pressure drops in 
tests 1 and 2 were 8.6 and 7.2 lb/ft”, Respectively. UEetecione.. tie 
seal design can be expected to influence the response. All further 
epemimental tests use the thick seals, 

Another important observation, which is not in the simulation, is 
the standing wave that is observed being created by the seals. This has 
the effect of causing local fluctuations in cushion volume and the exit air 


gap other than the normal rise and fall of the water level. 


3.0 Comparison of Simulation and Experimental Transients 

The transient response curve for the linearized equations is com- 
pared to the test data in Fig. 3.4. The natural frequency for the linear- 
ized equation of w= 9,2 rad/sec is the same as in Model A. The 
curves of lower damping are representative of the real situation where 
the faired seal (test 2) is used. This has a damping ratio on the order 
of &= 0.3. The experimental response (test 1) has a = 0.4 for the 
first quarter cycle of the nonlinear response. The addition of non- 
ImiMearity im tne exit gap area causes more cr enoot with the same value 
of €= 0.3, than was found in the totally linearized model equations. The 


woinme ot ¢ = 0.3 is carried through to the frequency response analysis. 
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Another result of the Model B is that the addition of the air mass flow 
has little effect upon the natural frequency and critical damping when 
compared to Model A without air flow. This is a very likely result 
Since the time constant for filling the air cushion region is long com- 
pared to the wave frequency. 

The last comparison to be made is the acceleration and pressure 
responses, Fig. 3.5, of the digital simulation and experimental tests. 
Agreement appears to be quite good except for the higher frequency 
oscillations which are thought to come from the bubbling action of the air 
through the water around the seals. Also some impact forces may be 
caused by the water slapping the seals. These localized dynamic 


responses are not accounted for in the Model B. 
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Chive isi 4 


SES FREQUENCY RESPONSE TO REGULAR WAVES 


fa. Introduction 

With the conclusions of the transient response the following chapter 
presents the frequency response analysis. Weave motion is the principal 
dynamic response to be analysed, Since it can be easily compared to 
experimental data. Information on acceleration and cushion pressure is 
also presented. The behavior in regular waves was measured in the tow 
tank with the SES craft constrained to heave only and operating at zero 
forward velocity. The simulation model is presented first and then com- 


Pee to time Experimental results. 


4,2 Model C Frequency Response 

The modeling of the heave dynamic frequency response is primarily 
based on linear theory main certain nonlinear effects also included. The 
linearized equations of Model C, as developed in Chapter 2, are included 
in (19). The third order system equations are the same as given in 
Model B except for the following additions: 

1. The external disturbances considered are the wave pumping (17) 
and the wave elevation at the bow and stern seals, (15) and (16). 

2. The principal nonlinearity that is included in the systems 
equations is the gap exit area (9), (10) and (11), wherein the leakage area 
is eno wied to vary as a function of the craft seal elevation and wave 
height. The model considers a nonflexible seal attached to the bow and 
stern. With the craft stationary in the x-y plane, air leakage exits from 
both seals. When the craft is moving forward then the principal leakage 


is from the stern seal only. 
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3. The body suspension formes acting on the craft, aside from the 
pressure force, represent the hydrostatic and hydrodynamic effects due 
to body-fluid interaction (including wave-induced effects) for the sidewalls. 
Estimates for added mass and damping come from the transient analysis 


emia [5]. 


4.3  Heave Suspension Performance 

In order to illustrate the nature of the craft responses to waves, 
digital computer simulations and experimental time histories are 
presented. The results of typical Simulations appear in Fig. 4.1. The 
responses to regular waves illustrate the magnitude and phase of the 
“heave, acceleration and pressure to the input. The constraint on the 
leakage area given by AY = 0 1S presented to illustrate the basic nature of 
the outputs. All motions are quite regular and the variation in gap area 
fmerie only nonlinearity. The behavior of the vertical CG acceleration is 
approximately peeoortional to but slightly lagging the pressure as the 
wavelength increases. The heave mean position, though not shown in 
Fig. 4.1, shifts downward as the gap area increases until a new equilibrium 
position is attained. 

Figure 4.2 illustrates typical time histories of the motion found 
during the experiments in the tow tank. Several comparisons between 
the simulations, Fig. 4.1, and experiments, Fig. 4.2, are made: 

1. Note the nonlinearity in the experimental heave response. The 

air leakage at the bow and stern appears to be ates the majority 

of the nonlinearity. The simulation model of the air gap leakage 
ssemmore linear. 

2. The experimental heave motion shows significant flattening of 


the response near the peak, which does not appear in the linear 





| Se a 
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Fig. 4.2 Experimental Heave Motion Time Histories fOr Sion 





Simulations. 

3. The heave motion lags the wave input and is fairly consistent 

berveen Experiment and simulation. 

The pressure and acceleration time histories from several experi- 
mental results are presented in Fig. 4.3. The high frequency oscilla- 
tions in the pressure and accelerations are felt to be caused by the air 
escaping under the seal due to the bubbling action which transmits pres- 
Pure pulses. likewise, they appear in the acceleration. The high 
frequency oscillations occur when the craft is reaching its peak in heave 
and the air exit gap opens. These high frequency oscillations, being due 
. to local effects around the seals, should be given careful attention in 
future studies. The high frequency oscillations, though, appear to have 
little effect in the overall heave motion prediction when comparisons 
between simulation and experiment are made on the macro level. In 
addition, the solid seals attached to the craft are transmitting the forces 
due to the wave impact on the seals. A flexible seal should result in 
Fower wave impact force. 

In order to give a better picture of what is taking place overall, let 
us analyse the heave dynamics over a range of wavelengths. Lhe Menace 
motion (rms), Fig. 4.4, from the digital simulation indicates a behavior 
somewhat similar to that indicated by the limit value given in (23), with 
near zero response at values of 


7X» 


oO” fwd 
jibes! IZ | = TT sin s (23) 
eee = |. 0. = = and =. Therefore, the minimum heave response will 


occur when the wave pumping is also a minimum, and likewise the maxi- 


es 2 ; 
mum response occurs close to A/éz = Zs 3» =e and 7. The small differ- 
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ence between the simulation and the values in (23) is due to the effect of 
the natural heave frequency which is near A/£= 0. 38. 

The large response characteristics at short wavelengths, which 
are due to the direct influence of wave pumping as well as to the reson- 
ance effect, result in relatively large heave acceleration (rms) values as 
indicated in simulation results, Fig. 4.5. This condition, which is 
generally known to occur in air cushion vehicles, is often referred to as 
"cobblestone" oscillation. It is evident when considering the operation 
of these craft over short wavelengths. This effect, though, was not 
verified from the experimental data, but is recommended for future 
verification. 

The results of the heave motion and acceleration from simulations, 
Pigs. 4.4 and 4.5, agree with the simulation studies in [1]. 

With the basic results of the heave motion presented, the simula- 
tions are now compared to the experimental test results. The heave 
mMe@tion frequency response data from several series of experimental 
tests are summarized in Fig. 4.6. Data for the response to waves of 
Henctinn = 1.5 - 12 ft, and height ae 2 in, indicate that the model 
exhibits the same general response as the experimental results. This 


is the minimums and maximums along with the increasing heave 


) 
response amplitude with wavelength are similar. Experimental data 
for wavelengths 3.1 - 4.6 were not available, Fig. 4.6; therefore, the 
results are not conclusive in this he gion- 

The most interesting result, Fig. 4.6, is the shift in the experi- 
mental heave response minimums as compared to the model. The two 


“primary minimums are at A/£ = 0.35 and 1.47. At the lower wave- 


lengths another minimum appears near A/& = 0.22. 
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In an effort to account for the differences between experimental 
and model results, several ideas were put forward based upon observa- 
tions of the cushion waves. With the lack of data for A = 3.1 - 4.6, the 
first assumption concerning the response in this area was to assume the 
heave response followed a smooth curve. This would then indicate that 
the primary minimum heave response shifts to the left and right of the 
model values A = 0.5 and 1.0. 

In order to determine what may cause this shift, obServations of 
the wave profile photographically indicated some differences between the 
external and internal waves. Further experimental testing was per- 
formed by mounting wave probes internally and externally to the SES 
craft. Experimental results for various wavelengths, Fig. 4.7, indicate 
that the incoming wave is being altered in several ways. First the 
internal wave is lagging the external wave by 20-30 degrees at the lower 
wavelengths. Secondly, secondary waves are appearing. Thirdly, the 
internal wave is amplified at data values A = 3.12', 6.5', and 9', and 
Slightly attenuated at the other wavelengths. This indicates that the 
seals and, to a lesser degree, the sidewalls are creating an internal 
Standing wave. For this reason, the minimum heave response would be 
expected to shift. The amount of shift would depend upon whether the 
standing wave or the external wave is dominant. The standing wave, 


having its own characteristic eigenvalue, could cause a shift in the 


system's natural frequency response away from the original model values. 


There is some speculation that a dip in the heave response might 
occur in the region A = 3.1 - 4.6 where experimental data were not avail- 
able, Fig. 4.6. If this were true, then the results would indicate that 


the minimums have shifted to the right of model values. Again, the 
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internal Standing wave would be the main contributor to such a shift. 

Since no further model simulations were performed that could 
account for the wavelength shift and experimental verification of the 
heave response at the missing wavelengths, this is left for future inves- 
tigation. Additionally, investigation with the craft moving at a forward 
velocity is needed. Sufficient questions have arisen that will hopefully 
induce others to further the research. 

Additionally, data comparisons of heave response to wave heights 
of 1, 2 and 3 inches, Fig. 4.8, indicate that heave damping increases 
with higher wave heights for A/2< 1.0. The maximum and minimum 
heave response to waves of varying height appears to shift toward the 
model values as the wave height increases, Fig. 4.8. This might be 
expected in that the external wave dominates any internal standing wave. 

This concludes the experimental tests conducted and the extent of 


verification of the simulation. 
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SUMMARY AND RECOMMENDATIONS 


5.1 Summary 

Theoretical model simulations and experimental craft tests for the 
See degree oi freedom heave motion ofga surface effect ship have been 
conducted. The theoretical third order linearized model dynamical 
equations consider the following: 

1. A constant air mass flow into the cushion. 

2, Quadratic exit gap orifice restriction for air flow out of 

the cushion, 
3, Conservation of momentum of the SES rigid body in the 
Mericaleanade. 

4, Conservation of air mass inside the cushion. 

‘Oo, Polytropie relationship between pressure, air mass and 

air volume. 

6, Stationary model in x-y plane, 

7. External force inputs and regular wave disturbances are 

the principal transients considered in analysis. 

In order to simplify analysis and to verify the system's responses, 
linear analysis has been developed and the basic third order characteristic 
equation derived, A study of the linear model characteristics under 
external force inputs has been conducted for the SES (4/b = 6.5) with 
constant air flow. In the best fit of the transient response simulation of 
the model and the experimental data, the natural frequency and critical 
damping ratio of 9.2 rad/sec and 0.3, respectively, were determined. 
From theory based upon Stoker's formulation [5], the natural frequency 


and damping, based upon an oscillating free surface pressure field, were 
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also found to be in good agreement with the experimental results. The 
transient response characteristics of Model B (with air flow) did not 
appreciably alter the response of the craft when compared to Model A 
(without air flow). This was due to the time constant for filling the 
cushion being an order of magnitude greater than the system's natural 
heave frequency. The experimental results verified this. The value of 
damping was determined to be the critical measure in estimating the 
amplitude of the transient response. As would be expected, the increase 
in damping helps to decrease the effects of external disturbances. 

The outcome, due to variations in seal design, was noted where the 
effective exit gap length under the seal influences the airloss. The 
solid seal design maximizes the effects of the nonlinear characteristics 
of the exit gap area, which, in turn, influences the heave motion, 
acceleration and cushion pressure. These results would likely be 
attenuated with flexible seals. 

The sensitivity of vehicle performance to wavelength and wave 
height were studied with the following primary results: 

1. Prediction of heave amplitudes from theory compared favorably 
ieee eperiiirchital results. The relative amplitude to wave height 
ratio was found to decrease with increasing wave height for A/£ $s 1.5 in 
PMiesexperiinental tests. This is indicative that damping increases as the 
wave height increases, 

2. The predictions of minimum wave pumping from theory, 
assuming the Froude-Krylov hypothesis, should occur near A/2£= 1, oe 
and Pe Experimental evidence, though, indicates that at low speeds the 
primary wave is aitered by the seals and sidewalls. 


3. AS aresult of seal and sidewall effects upon the waves under the 


06 
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cushion, the points of minimum and maximum heave motion are altered. 
Experimental evidence shows that the minimum wave pumping is shifted 
to the left and right of the two primary theoretical values. For this 
particular test craft, the minimums occurred at A/£= 0. 38 and 1.47. 
Comparisons of the inside and outside wavc profiles from experimental 
tests indicate a general lag in the wave under the craft whcn compared 
to the outside wave. Also the amplitude of the waves under the cushion 
is attenuated in some wavelengths while others are amplified at other 
wavelengths. These effccts are not considered in the simulation model 
and therefore account for the major difference between theory and experi- 
ment, 

4, The sensitivity of acceleration levels to external disturbances 
was studied in the simulation model. The accelerations are 180 degrees 
out of phase to the heave motion. The acceleration (rms) response over 
a range of wavelengths in the simulation gives the ''cobblestone'’ appear- 
ance with relative magnitudes increasing at shorter wavelengths. No 


experimental tests were conducted to verify this result. 


Turning now to the experimental results, high frequency oscilla- 
tions of the cushion pressure and accelerations were observed, This 
phenomenon occurred when the exit gap area was opening and closing. 
The air appcars to escape under the seals with bubbling action. This 
effect is not apparent in the simulation models since a smooth opening 
and closing of the exit area is considered. These high frequency 
oscillations, though, appear not to alter the basic craft hcave response. 
Also, there may bc somc water impact forces on the solid seals which 
could account for some of the high frequency accelerations, 


In conclusion, the third order linear model equations are fair 
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representations of the model dynamics of the craft and air cushion. The 
principal differences in heave response between theory and experiment 
would appear to be attributable to the water dynamics under the cushion. 
The next step in the modeling must allow for this seal-water interaction 


and the resultant waves under the cushion. 


9.2 Recommendations 

1. A detailed study of the waves inside the cushion along with their 
cause, i.e., Seals and sidewalls, and their effect upon the cushion 
dynamics is needed. 

2. Modeling techniques as a result of 1 could, therefore, predict 
the system dynamics to a greater degree of accuracy. 

3. Regions of applicability of the Froude-Krylov hypothesis should 
be established, 

4, Although the quadratic orifice equations appear to predict fairly 
accurately the air mass flow from the cushion, measurements of the time 
lags induced by the air flow around the seals and their effects upon the 
overall dynamics should be considered in a more detailed study. 

9. A study of the contribution to damping and added mass by the 
sidewalls would be a valid consideration not studied in this thesis. 
Slender body theory is suggested as one approach. 

6. Further experimental model tests should be conducted with the 
craft moving with a forward velocity. 

7. Acceleration amplitudes between theory and experiment need 
verification. | 

8. After completing the wase study, the pitch motion of the SES 


may be investigated. 
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APPENDIX A 


EXPERIMENTAL FACILITIES AND TESTS 


A.1 Introduction 

The Sio@vithelenpth to beam ratio of 6.9 was Selected for study in 
the experimental program of this thesis. Since little open and detailed 
literature is available on the motion (heave) characteristics of this 
particular SES, a test program was formulated and carried out. Results 
from the Static, transient, and frequency responses were used.in the 


Simulation study. 


m2 Test Program 


For the purpose of this thesis, the experimental test program was 
designed to accomplish measurements of (1) the static stiffness coefficient 
(displacement/load ratio) of SES, (2) the transient response to a step input 


in force, and (3) the frequency response of the SES to a series of regular 


Waves, 


=. lest @memmeamderacility 

Figure A.1 is a photograph of the SES craft in the tow tank with the 
tow carriage and associated instrumentation, The basic test craft 
dimensions and equilibrium testing conditions are given in Table A.1. 
The test craft's sidewall and deck were constructed from 4-inch plexi- 
glass, The sidewalls and seals at the bow and stern were shaped (faired) 
from styrofoam (Fig. A.2) and attached to the feet eaeutees Note: the 
seals are inflexible and positioned for a particular equilibrium state. 
One advantage of building tne test craft from plexiglass is the visual 


Observation of the free surface inside the cushion. Photographs of the 


wave profile were taken from the side, 





Fig, A.1 SES Tow Tank Model 
Attached to Carriage. 





Fig, A.2 SES Bow and Seal Configuration. 
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TABUE A. 1 


Seo MODE Le Gran ACT Rintsl lcs 


Cu Shvonmibe not hye ity) ue. 5... ues 6.26 

CUES, 0 ON MB Toye icc) 0) CR) Nees aaa ear ee nner erry ar 0.96 

Mepeth BedmehealiO . 26... ee ee eee 6.9 

Memopie@veradlis Lit): 2. .ccle saris wale. 6 eee ee 6. 66 

ipeumim@nmeneai I(it) 2... 36605 cs ca Bla ehes es 1.19 

Cushion Area, A (ft”) Siete eM MC elec wa tens ar oeeess 6, 00 

Sidewall Height, z, (it) ees. of eee ee 0,328 

Cushion Height, z, (ft) ek ee ee ee OFZ 

Sidewall Thickness, Y sy, (Ft) i) A, ate ee 0.115 

Vented! Center of Gravity st) fae. 2. ose ec 0.386 

Longitudinal Center of Gravity, (ft) ........ 3.27 from stern 

TO Ged Wiene ity UNO) scree, aay cts seners me ae «nae ee eo (2.0 

Specific Weight Density, (1b /ft®) seven enone : aw) 

Mie oupply PlOweinetes J... 26s 6 < eee Out te 3 min (Sa 
ee eee #2 19 ft, °/min ey 

[NE ISO? gem ne ec 2 oe 0.075 lb/ft? 

Water Density .....cccceeceeccecees —ae 62.4 lb/ft? 
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Pressurization of the cushion was accomplished by two pressurized 
Supply lines, which were 3/8-inch in diameter, located near the bow and 
Stern. Figure A. 3 illustrates the air supply and instrumentation 
involved, The air flow, hovering height, and craft weight may be 
adjusted independently. Nonequilibrium conditions were simulated by 
introducing regular waves, Several pressure transducers, wave probes, 
a heave rod, an accelerometer and an external flow meter were provided 
to accomplish quantitative studies, A chart recorder was available for 


recording each response. 


A.4 Static Response 


The static response for craft height versus load was conducted by 
varying the craft weight and measuring the vertical displacement from 
mve heave rod, See Fig, A.4for results. The observed stiffness was 
approximately linear over the operating range. The linearized spring 


coefficient (k, + ky) is equal to 435 + 10 lb/ft. 


A.5 Dynamic Response 

Transient responses were performed by removing weight from the 
craft and allowing it to reach a new equilibrium. This is equivalent to a 
Step input inforce. Transient performance data for heave, acceleration, 
and pressure were obtained. 

The frequency response testS were performed by pasSing waves of 
various lengths and heights past the stationary SES. Response data were 
recorded for heave, acceleration, pressure, and the wave profile inside 
and outside of the cushion. Due to the noticeable difference between the 
free surface inSide and outside the cushion, two different means of 


measuring this reSponse were devised. Initially, photographs were 
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With Air Flow 


No Air Flow 


Weight 
(1b) 50 


k +k, Sosa one) i 


40 


0 Oso 


Vertical Displacement z, (in) 


Fig. A.4 Experimental Stiffness Coefficient of SES. 
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taken of the inside and outside waves through a viewing window from the 
Side of the tow tank. From these photos it 1s possible to observe the 
variance in the waves, but quantitative data analysis requires better 
photography. The second means of obtaining the time history of the 
cushion wave profile was by attaching wave probes at midship to the 
Sidewall, both inside and outside. The wave profile, relative to the 


craft, was obtained, 


Discussion of the results of these tests is found in the main text, 


65 





P92) 27) 10 BY Dee : 


oot Wino AD DAMPING FOR OSCILLATING FREE SURFACE 


Figures B.1 and B.2 were used to estimate the cushion added mass 


and damping. 


Added Mass Coefficient from [5]. 


Wave Number( oe) 


Pasi. 1 
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APPENDIX C - COMPUTER PROGRAM FOR MODEL C 


HEAVE DYNAMYCS OF SURFACE EFFECT SHIF 
DIGITAL SIMYLATION USING DYSYS (SUBROUTINE) 


VAR TARLES LISTED AS FOLLOWS 
A= ADnED MacS L8 SEC2,-T 
ABr CUSHION ARCA (FT«4>) 
AEFEQUYL + AIO FXIT AREA FROM SEALS (FT«e?) 


AEB EXT AREA aT BOW (FT #e2, 
AES= EXry7 AREA aT STERN (FTee2?) 


= DAMPING (; E*®SEC/FT) 
COG TOS mE NTRAICE AND FXYIT t ESISTANCE CCEFe (O.6) 
FR= AIp TEMP (¢) 
G= GRavITY 22,2 FI/SECw*? 
Me= EQHILIBSrU4 HEIGHT OF CUSHION (FT) 
{= 
A 


= AMAJENT ain PRESSURE ABS (LB/FT2) 
Me-BeTIGelON PSESSUFE ARS(ILR/FT2) 

PK= POLYTRGOFIC CONSTANT (1,2) 

GAs MAss AIR Fy OX RATE INTO CUSHION (SLUG/SEC) 
Rw= WATER DENS TTY (LB*SECHROsFT ome ) 


RwWS= AIR DENSI Ty (LB*SEC#*2/F Taw) 


TN= WAVE HEIGHY (FT) 
= FORWARD SrPeEO (FT/SEC) 
w= CRaFT WEYGuT (LB) 
WE FREQ OF WAVES (RaDyGeC) 
*L= WAVE LENGTHY (FT) 
KM= CRAFT MASS ;(L5eSEC2/FT) 
¥W= ENCOUNTER FRED FOR wEAD SEAS (RAD/SFC) 
¥B= CUSHION LENGTH CFT) 
YB= CUSHION wiInTr (FT) 


YSW= SInEWALL wyOTH (FT) 


Cn= SEAL PHETGHY (FT) 


ZIN= AVE, wAVE WEIGHT UNDER CUSHION (FT) 
ZOB= WAVE HEIagny AT BOW (FT) 
ZOS= WAVE HEIGHT AT STERN (FT) 


MS 2aePOCSITInN OF CRaFT 

Ve) Sees VorticCAL VELOCITVy OF CRAFT (FT/SEC) 

Y(3)= aIR MASS IN CUSHION (SLUG) 

Y(S)=Pc“xCuSHyON PRESSURE (LRB/FT2) 

F(2@)" aCCELERATION OF CRAFT (FT/SECHa?2) 

F( 3) a~IR “MASS FLOW RATE (SLUG/SEC) 

¥(9)= yNPUT wAyE HEIGHT (FT) 

Y(i4)® GAP HEIGHT FOR BOW ALD STERN (AVE) (FT) 
SUBROUTINE ¢QSIM 
COMMGN Ty dT, YCA2)sF (2OdsSTIMESFYTIME,MEWDT,IFLRT YN 
IF (nEWOT) 42352 ; 


CONSTANT FeR EQUILTSR Tuy GF EXIT AIF GAP UNLER SEAL (lesa) 
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REArn( Rei) YBsYBsYSW,ZSW We F Re GAsHlgbhtfs As Bs TR aL yUs PQs AB y PMO 


FORMAT (RE 1H ,3) 
FORMAT(LE4¢ 2) 
RWEE2 64 

TNEYIN/S (12.0% un2°O) 
WMawys32e2 

G*32.2 

ABzypeYB 
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PAt3416°R 

TKeS3¢34u(4-0eA+FR) 

RWSepAaysTK 

WMAsWhli+A 

WFR (DeM*3,1,E*G/HL) OOS 

WHEWF +(WRaeoe@) eU/G 

CO=ea,6 

CE=A.6 

PHASEK3¢ 142 eXB/hL 

RETIN YReWLE GIN (PHASE )/(3,142*AB) 
F(A yp sV¥ eo DAAC AS 

PIiga, VCS {7 wha 

F dete <7 

PK=4q.? 

PTaHQeRWs(PyeePO)4109 

PKR3aCOK (Pe Fa RAS /G) eA SHAE Ge 5/ (POPPA) eeBeS 
PRK4=xCOe(P,PaRYS/SG) mare Sane GeYBue (POPPA) 4005 
ClleaResPT 

CiSsH? *RwWshP aC 

DIie-sety/PT 
A@lase(RWeKX By YSWH? 0 O4nB¥C{ 4) /WMA 
Aetee S/W i 

AC3nceAGECITS/HMA 

GElajecd/WweaA 

Beea(RA*XB*ySnee oVeaReDIie) SWMA 
A3BLanPK3eC1 ye PKG (m1 6 S4C 11 /RW) 
A338 ae" K3aCi ger KaeC13/Rb 
B3ene?K3ePloerKan (Le atDiesRw) 

WAHT awe T 

ZOBSTNxSIN( wwTePRASE ) 

ZOSSTN«SIN( wWTePHASE } 
ZINSRYSIN(WwT) 
Y¥(SpaCii¥y(4)+C13*Y(2)4+012*ZIN 
YODya=TNeSIy (wl) se. 4 
YCleyp="¥ (4) /B 9°40 

Y¥C1l4q)="¥(2) ,6°52 

YOlaypeeF (PP) 710.8 

YOL3Z;SYC3) Syed 
YCL4)SAES(2C  O* YB) *12.40 
Yllgysv(5)/40°0 

AEBsyix(Hy+ (PO ePA)/RW470B¥Y(4) ) 
AESsyhxlhie( PU =PA)/RW+70SHY(41)) 
IFC aARReLE,t 6) AEBER LD 
IF(CARSeLEeY od) AESEA,S 

ACR ALR+AES 

PK3=aCOMR( A eV URWSSG) MEMOS e AE ee S/( POPPA) wae S 
ASlaeFPK3aC1L ePK4e (me O4+C1 L/RW) 
A3S3aeF K3¥C1iaePK4*C13/Rh 
KSGae"K3eN oer Kha (Le +D1LP/ Rw) 
FliypxYle) 

FCC ychliny ly) +A2C#Y (2) FAC 3*Y (3) 4R2O1eF 1 +Be2eZIN 
F(3)aABL¥Y (4) tAR3™Y( 2) 482208 ZIN 
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